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Parameter estimation is of fundamental importance in areas from atomic spectroscopy and atomic 
clocks to gravitational wave-detection. Entangled probes provide a significant precision gain over 
classical strategies in the absence of noise. However, recent results seem to indicate that any small 
amount of realistic noise restricts the advantage of quantum strategies to an improvement by at 
most a multiplicative constant. Here we identify a relevant scenario in which one can overcome this 
restriction and attain super-classical precision scaling even in the presence of uncorrelated noise. We 
show that the quantum improvement can be significantly enlarged when the noise is concentrated 
along some spatial direction, while the Hamiltonian governing the evolution which depends on the 
parameter to be estimated can be engineered to point along a different direction. In the case of 
perpendicular orientation, we identify a state which achieves optimal asymptotic precision scaling. 



I. INTRODUCTION 

Estimation of an unknown parameter is essential across 
disciplines from atomic spectroscopy and clocks [1-3 to 
gravitational wave-detection [4 . It is typically achieved 
by letting a probe, e.g. light, interact with the system un- 
der investigation, picking up information about the de- 
sired parameter. As seen in Fig. [l] a metrology protocol 
can be understood in four main steps [5^, ^6] : i) prepa- 
ration of the probe, ii) interaction with the system, iii) 
readout of the probe, and iv) construction of an estimate 
of the unknown parameter from the results. Steps (i)-(iii) 
may be repeated many times before the final construction 
of the estimate. 
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FIG. 1. General metrology protocol where a known probe 
state evolves according to a physical evolution depending on 
an unknown parameter uo. After sufficient amount of data is 
collected an estimate for the parameter is constructed. 



The estimate uncertainty will depend on the avail- 
able resources, here the probe size N and the total time 
T available for the experiment (other choices are possi- 
ble [7]). By the central limit theorem, for N uncorrelated 
particles, the best uncertainty scales as l/y/vN ^ where 
V = T/t is the number of evolve- and- measure rounds. 
This bound is known as the shot-noise or standard quan- 
tum limit (SQL). By making use of quantum phenom- 
ena, a metrology protocol may surpass the SQL, reaching 
instead the limits imposed by the quantum uncertainty 
relations. For probes of non-interacting particles, the 



best possible scaling compatible with these relations is 
l/(^JvN)^ known as the Heisenberg limit. 

Without noise, the Heisenberg limit can be attained 
using entangled input states, e.g. Greenberger-Horne- 
Zeilinger (GHZ) states for atomic spectroscopy [8 . In the 
presence of noise however, the picture is much less clear, 
as the optimal strategy depends strongly on the model of 
decoherence considered. Nevertheless, the SQL has been 
significantly surpassed in experiments of optical magne- 
tometry O [10] , which proved that some sources of noise 
can be effectively counterbalanced [TTJ|12^. However, un- 
less one can keep improving the interaction strength or 
readout efficiency with probe size, e.g. increasing the op- 
tical depth with the atom number in atomic vapours, 
destructive effects of uncorrelated noise are bound to 
dominate at high N . In this regime, a number of no- 
go results exist, demonstrating that for most types of 
noise acting independently on each probe particle, an 
infinitesimally small amount of decoherence limits any 
quantum improvement over the SQL to at most a con- 
stant factor p^l fTSHTS] . In particular, any noisy evolution 
described by full rank channels p!5HT8 belongs to this 
class. This is arguably the most likely kind in experi- 
ments, suggesting that the precision is always bound to 
scale classically for large enough N . 

Here, for a frequency estimation task undergoing un- 
correlated noise, we show that the restriction to SQL- 
like scaling can be surpassed even when the process 
is particle- number independent. This is achieved by 
optimising the duration t of the evolve- and- measure 
rounds, for a configuration where the Hamiltonian and 
the noise have preferred spatial directions transversal to 
each other. The corresponding channel is full rank, but 
we find asymptotic precision scaling of To our 

knowledge, this is the first time the SQL is surpassed for 
full-rank noise. We demonstrate the scaling to be opti- 
mal by constructing an upper bound on the attainable 
precision and identifying a state which saturates it. We 
also analyse deviations from perfectly directional noise. 
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The asymptotic scaling is then again restricted to SQL- 
hke. However, for smah deviations, we observe a much 
higher precision gain than for perfectly parallel noise. 



II. MODEL FOR NOISY FREQUENCY 
ESTIMATION 



We consider a Hamiltonian H = f Xl^i^s^ where 
is a Pauli operator acting on the /c'th spin- 1/2 particle 
(qubit), and uj is an unknown frequency to be estimated. 
To account for noise, we model the evolution by a master 
equation of Lindblad form 



dp it) 
dt 



n{p) + c{p). 



(1) 



Here, H (p) = —i [H, p] describes unitary evolution and 
the Liouvillian C{p) describes noise. We consider uncor- 
related noise, such that C = "^j.C^ and for a single qubit 
we have 



[P- 



aya2P(T2 



a.alpa^] 



(2) 



where 7 is the overall noise strength and (^x,y,z ^ with 
o^x-^<^y-\-az = 1. For az = l this describes the situation 
considered by Huelga et al P, namely dephasing along 
the direction of the unitary, while ax = l corresponds to 
dephasing transversal to the unitary. For 
1/3 we have white noise, i.e. a depolarizing channel. 

The uncertainty 5uj in the estimate of uj can be ex- 
pressed in terms of the quantum Fisher information [19] 
(QFI) J^{puj) of the probe state p^j after evolution. Ac- 
cording to the quantum Cramer- Rao inequality [20] , 



6ujVf> 



(3) 



a bound achievable asymptotically ioi v ^ 1 [STJ [22] . 
One can explicitly solve the master equation ([T]) , obtain- 
ing a map for the corresponding channel (see Section IV 
for details). Applying the map to a given input, one 
obtains p^^ and the QFI is computed through the diag- 
onalization of this state. Alternatively, bounds on the 
QFI can be computed from the Kraus representation of 
the channel [El [15]. We write = J^{puj) for probes of 
particles. 

For a general input state, it is a difficult task to com- 
pute the QFI, as the size of p^ grows exponentially with 
A^. However, for a GHZ state, 



\GHZ) 



V2 



|1)^ 



(4) 



the calculation simplifies dramatically (see Section |IV[ ) , 
and we are then able to optimise J-'N/t (and hence Su) 
for large N. To obtain results valid for general inputs, 
we resort to bounds on the precision. A first indica- 
tion that performance is better under transversal than 




FIG. 2. Cut through the space of all valid qubit maps. The 
top point corresponds to the identity channel X[g] = and the 
disk circumference to (extremal) unitary channels Uujt- One 
can think of phase ojt as increasing along the angular direc- 
tion and decoherence jt along the radial direction. The spi- 
rals represent trajectories of channels for parallel S^{t) and 
transversal S^{t) noise, approaching the completely mixing 
channel at the disk centre. Arrows are separated by fixed 
time steps. We see that S^{t) loses coherence slower, being 
always closer to the boundary, and also has larger angular 
'speed' \d8^/doj\. The classical simulation method l^ pro- 
vides a geometric bound on estimation precision based on the 
distance to the boundary: SljliVT>^ {e{tYt)/N . For asymp- 
totic A/", the bound is determined by the behaviour for small 
t at the start of the spirals. We get ez{tY — 2'^ /t + 0{l) and 
e.(t)' = (7't')/12 + 0(t^). 



parallel noise is given by the geometric Classical Simula- 
tion method [T^ (see Fig. [2|. We obtain tighter bounds 
by adapting the finite- channel extension (CE) method 
of Kolodynski et al. [23] to allow for optimization over 
t (see Section IV). As for previous asymptotic meth- 
ods [m [15] without t-optimization, we cannot apriori 
guarantee our bound to be saturable. However, for par- 
allel noise (a^ = 1), our bound is known to be tight in 
the asymptotic N limit both with [24 and without [23] 
t-optimisation. As we find below, the bound is also tight 
for transversal noise. 



III. NOISY METROLOGY BEYOND THE SQL 

Parallel noise. To understand the role of noise we com- 
pare parallel and perpendicular noise, as well as noise 
which is directional but not fully concentrated on any of 
the axes. We start by parallel dephasing, which has been 
studied before [T. 

For a single qubit, the QFI for an optimal input state 
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FIG. 3. (a) Precision scaling for bounds and GHZ strategy. The bound for parallel noise with optimisation of the single-round 
duration is shown (dotted), as well as the bounds for transversal noise without (dashed) and with (solid) optimisation. Without 
optimisation (taking t optimal for a single qubit) the asymptotic scaling is still SQL-like, but with optimisation super-classical 
scaling is maintained. In the latter case, the bound is saturated by the GHZ-state strategy (dots). The thin line shows 
scaling for reference, and the borders of the upper and lower shaded regions show SQL-like and Heisenberg scaling respectively, 
(b) Bound (solid) and GHZ-state strategy (dashed) for noise with a dominant transversal and a small parallel component, 
OixjoLz — 99. The scaling is super-classical for moderate N ^ while the 2;-component determines the asymptotic scaling which is 
hence SQL- like. The CE bound approaches (lower thin dotted) and the GHZ-state (middle thin dotted). 

The finite- A/" CE bound for parallel noise with the same strength is shown for reference (upper thin dotted) as well as for 
transversal noise (thin solid). 



is given by = e~^^^t^. In the noiseless case (7 = 

0), as expected, the longer the probe state is allowed to 
evolve, the more information can be extracted about the 
parameter. In the noisy case, after a certain time the 
dephasing process wins over the unitary evolution and 
the extract able information is degraded. For a classical 
strategy using N independent qubits = NJ^"^^. The 
best performance is obtained by optimising the evolution 
time to minimize 6uj. For a single qubit t^^^ = I/27, and 
the SQL for parallel noise is then 



5uj VT > 



2je 



(5) 



It has been proven [14 that when quantum strategies 
are allowed, the precision is instead bounded by 



SujVf> 



1 + (e27^ - l)N 



> 



N 



(6) 



where (7, t) = (e^^^ — l)/t. The tighter bound coincides 
with the finite- CE method [23 for this channel and 
asymptotically reduces to the weaker bound (also found 
by the method [15 of Fig. Optimising over t we find 
t°P* =^[Ar]/27, where w[N] = 1 + and W[z] is 

the Lambert W function. Asymtotically, t^^^ approaches 
zero as Hence the asymptotic scaling of ([g) is 

dictated by the small-time behaviour of 0^(7, t), which is 
fully determined by the channel geometry (see Fig. [2): 



,(7,t) = 6,(^)^^ = 27 + 272^ + 0(^2). 



(7) 



From (|6| we see that the asymptotic scaling is SQL-like. 
Comparing to ([5|, it becomes clear that quantum strate- 
gies provide at most a constant factor improvement of 
y^, as found earlier [1] |14] . The stronger bound of (|6| is 
shown in Fig. [sj^a). Asymptotic SQL scaling is evident. 



The asymptotic SQL-like scaling is known to be sat- 
urable with spin squeezed states [23l [24] . For compari- 
son with transversal noise, we note that for a GHZ input 
state, the resulting bound on precision is identical to the 
SQL ([5|. A GHZ strategy is thus useless for parallel 
noise. 

Transversal noise. We now turn our attention to per- 
fectly transversal noise. While we do not have an analyti- 
cal expression, we can efficiently compute the finite- CE 
bound and determine the optimal evolve-and-measure 
duration for each N (see Fig. [sja)). 

Without optimisation, the asymptotic scaling is still 
restricted to SQL-like, as indeed it must be since 
transversal noise corresponds to a full rank channel [15]- 
18 . However, the channel geometry allows for an ar- 
bitrarily large asymptotic quantum improvement (see 
Fig. [2|. At short times, the asymptotic SQL constant 
factor Ca;(7,t) obtained by classical simulation is 



^2 



(8) 



As opposed to ([7|), Cxi'y^t) ^0, as t^O. Hence, unlike 
for parallel noise, the asymptotic quantum improvement 
factor is geometrically unbounded. 

Optimising the finite- CE bound over t for each A/", 
we find that super-classical scaling is maintained (veri- 
fied numerically up to A^ = 10^ for 7/cj between 0.001 
and 100) and follows an asymptotic behaviour very well 
described by 



5ujVf> 



Ar5/3 • 



(9) 



For large A/", 0, since otherwise the restriction to 

SQL-like scaling applies. From our numerics we obtain 
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that t°P* = (3/7Ar)V3 and 



o2/3 

2 



(10) 



Thus, at the level of bounds, the t-optimised quantum 
strategies provide a scaling rather than a constant factor 
improvement over classical schemes for transversal noise. 
To confirm this scaling, as we have not yet demonstrated 
that the bound is tight, we examine a specific strategy 
based on a GHZ state. 

We can analytically compute the QFI corresponding to 
a GHZ input (see Section IV). The expression becomes 
cumbersome with larger A/", but we can numerically de- 
termine and the minimum ^cj, for values of N up 
to several thousands. The result is shown in Fig. |3ja). 
Clearly, for the displayed values oi ^juj the GHZ state 
is very nearly optimal. What is more, the GHZ strategy 
shows no sign of returning to SQL-like scaling for large 

(verified for A^ up to 5000 and 7/cj between 0.001 and 
10). Based on this strong numerical evidence, we con- 
jecture that the finite- A" CE bound is indeed tight for 
sufficiently large A^, and that the asymptotic scaling is 
super-classical as predicted by ([9|. 

Intermediate noise. In a realistic implementation of 
a setup with transversal noise, most likely there will be 
deviations from perfect directionality, i.e. there will be a 
small amount of additional parallel noise. We can model 
this by taking ax = 1 — e, and = e. Once there 
is some z- noise, the asymptotic scaling must return to 
SQL-like. If this were not the case, then by starting with 
a bit of noise and adding x- noise until ax dominates, 
one could recover super-classical scaling. Hence, for a 
sufficiently large A", the precision would be improved by 
adding noise, which is clearly unphysical. Nevertheless, 
super-linear scaling can still persist over a large region of 
N. 

Following this reasoning, the asymptotic bound must 
be of the form 



5ujVf> 



.(7,e) 



A^ 



(11) 



with Cxzij^^) ^ Cz{ej,0) = 2€j. In fact, we expect that 
equality must hold, so that (11) is determined solely by 
the parallel component of the noise. This is because 
t°P*^0 for asymptotic A^, but for very short times non- 
commutativity effects seize to apply and our model is 
equivalent to a process in which unitary evolution and 
noise along each axis are applied sequentially in any or- 
der. Hence in the asymptotic limit, we can always think 
of perpendicular noise as being applied before parallel 
noise. As shown already, during the first stage scaling 
is super-classical, so that the asymptotic scaling is de- 
termined solely by the second. This is confirmed by the 
numerical results in Fig. [sj^b), where we see that the CE 
bound attains the asymptotic scaling y^zTe/iV. 

The same argument applies to the GHZ-state strat- 
egy. In this case, we can see explicitly what happens at 



short times. Expanding up to second order in t we have 
TN/t = NH - [{2N - 1)(1 - e) + m^e] Njt'^/2, and it 
is clear that for sufficiently large N the contribution from 
transversal noise is negligible in comparison with the 
parallel-noise contribution. The expression then reduces 
to the case of pure parallel noise with an effective noise 
strenght of 67 and asymptotically the precision scales as 
y^2eje/N. Thus, comparing to ([ll]), for e > the GHZ 
state no longer saturates the bound. In Fig. [sj^b) we 
explore the transition from super-classical to asymptotic 
SQL-like scaling for both the bound and the GHZ state. 
We note that before the transition to SQL-like scaling, 
GHZ appears to be nearly optimal. Since for small e this 
transition occurs at very high A", we conclude that for 
finite probes and small deviations from transversal noise, 
GHZ states offer precision enhancement very close to the 
maximal. Intuitively, we expect that the transition point 
Nxz{li e) to SQL-like behaviour increases smoothly to in- 
finity for e ^ 0. Indeed, from ([9|, ([To| and ( pT] ) we can 
estimate Nxz as the intersection of the super-classical and 
asymptotic SQL-like asymptotes, Nxzij^ e) ~ 8/(876^/^). 
Beyond A^a^(7, e), no significant gain is obtained by in- 
creasing the number of entangled particles in the probe, 
i.e. increasing the size of the entangled probe or using 
copies of disentangled probes of the same size leads to 
the same improvement in precision. 

Finally, we remark that while tends to zero for 
large A^, in practice the evolution time that can be im- 
plemented is bounded from below and can be set only up 
to a finite resolution t^^^. However, even with this limi- 
tation, large improvements over classical scaling may be 
obtained. In the case of transversal noise, one can opti- 
mise t for A^'s up to the point where =t"^^" and obtain 
the l/AT^/^ scaling of Fig. ^a) in this region. Beyond 
this point, one sets t = t^^^^ and the precision becomes 
restricted to SQL-like. Thus, the effect of t"^^" > for 
transversal noise, despite its different nature, has similar 
consequences to the effect of e > 0. Both lead to non- 
zero asymptotic constant factors, Cx{j^t) and €3^^(7,6) 
respectively, which can approach zero when t^^^ ^ or 
e ^ 0. More generally, when both e > and t^^^ > 0, 
the ultimate precision is limited by the bound (|6| with 
c.(7e,t™). 



IV. METHODS 

Solving the master equation. We write the map cor- 
responding to evolution under the master equation ([T]) 
for time t as a composite map of the form . Follow- 
ing Andersson et al. the single-qubit maps are then 
given by 



^u; (p) = ^Sjjdipdj, 



(12) 



where the di are normalised variants of the Pauli 
operators di = (T^/a/2, and all elements of the matrix S 
are zero, except ^oo = a + 6, Su = d -\- S22 = d — 
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5*33 = a — Sqs = 2C, S'o3 = —ic. For the application of 
known bounds for parameter estimation based, the map 
can easily be put on Kraus form [25] (see below). 
The coefficients a, 6, c, and / are in the Appendix [Aj 

GHZ input From the symmetry of \GHZ) and of 
the channel under permutation of the parties, one finds 
that the evolved state is diagonal in the GHZ ba- 
sis, that is the basis formed by all states of the form 
(|mi • • • mAr) ± |mi • • • mAr)) /a/2 with = 0,1 and 
\mi) = cri\mi). We can then parameterize each density 
matrix element (m^ . . . mAr|p0|m^ . . . m^) by the number 
m = ^ m^. It is possible to show that the diagonal terms 
are given by 



where A™(t) correspond to the optimal Xa/b in (16) ob- 



tained after performing the SDP minimization over h(t) 
for given t. After finding the optimal t = t'^^^{N) for each 
N in (17), we arrive at the final bound applied to all the 



t-optimised models considered in the paper, that is, re- 
placing t^t°P*(A^) in ^7^. As shown in Fig. (sja), the 
bound (17) can also be used to upper-limit precision of 



strategies not incorporating the t-optimisation. 

Crucially, as the second term in (17) depends on N 
through both A^^" and it is not true that in the 

limit N ^oo it is optimal to set A5 = as in the asymp- 
totic CE method [15]. For the purpose of this paper 
we have implemented a semi-definite program using the 
CVX package for Matlab [26]. 



1 



{m\p^\m) = -[d^a''-^^d 



(13) 



DISCUSSION 



while the anti-diagonal terms fulfill m' = = N — m^ 

and hence 



r{b-zcf-'''^f 



pN- 



'{b + icY 

(14) 

Each m-term has a degeneracy of (^). Using this 
parametrisation, the diagonalization of the 2^ x 2^ 
density matrix is reduced to the diagonalization of 
[A^/2 + Ij different 2x2 density matrices ([-J denotes 
the fioor function). 

Finite-N CE method. The channel extension method [151 
[23] uses the fact that allowing the channel to act trivially 
on an extended space can only increase the precision of 
estimation. Hence, for N independent channels 



max [£c 



(15) 

As a result of the extension of the input state space, an 
upper bound on the QFI can be obtained without any 
input state optimisation . The resulting minimisation 
problem can be cast as a semi-definite program [15] (see 
the Appendix Ib|. This is the finite-TV CE method [13] 



J^N{t) < 4min{7VA, + N{N - l)Xb}, 

h(t) 

subject to: A{t) > 0, B{t) > 0, 



(16) 



with Acj/5 real parameters, h(t) a rxr Hermitian matrix 
parameterizing locally equivalent Kraus representations 
of the channel, and A{t) and B(t) matrices depending on 
the Kraus representation. 

To accomodate for the optimisation over the evolve- 
and- measure round duration, we adapt (16) and max- 
imise over t the bound 



< 4 



N- 



'^™"-(^)+iV(iV-l)^ 



(17) 



Although recent results indicate that realistic noise 
prevents quantum metrology strategies from outperform- 
ing their classical counterparts by more than a constant 
factor, it is indeed possible to observe scaling beyond the 
standard quantum limit in the presence of uncorrelated 
noise. Adapting the classical simulation and finite-A/" 
channel extension methods [HI , we have shown that 
optimising the duration of evolve- and-measure rounds 
significantly enhances the precision. In particular, we 
have considered a unitary evolution with a well de- 
fined direction, and noise with a preferential direction 
transversal to the unitary. In this setting we have anal- 
ysed a frequency estimation protocol and showed that 
the GHZ state achieves maximal precision in the pres- 
ence of noise, providing a precision scaling of N~^^^. 
Furthermore, we have demonstrated that although the 
asymptotic scaling returns to SQL-like when the noise 
deviates from being perfectly transversal, the constant 
factor improvement of the precision can be significantly 
enhanced for small deviations. From a practical point 
of view, our results can be summarised in the following 
message: For a system described by the evolution model 
([1]), with freedom to choose the direction of the unitary 
evolution relative to the noise present, and for probes of 
many particles such that the optimal evolve-and-measure 
round duration is beyond the time resolution of the ex- 
periment, it is optimal to align the Hamiltonian with the 
smallest component of the noise and set the round dura- 
tion as small as possible. Doing so greatly amplifies the 
quantum enhancement of precision in noisy scenarios. 

We believe, our work opens an avenue towards useful 
quantum metrology protocols in realistic, noisy settings. 
To preserve a superlinear scaling in the noisy scenario 
there must be a well balanced trade-off between the size 
of the probe state and the time it evolves through the 
channel. We further hope to motivate novel work in this 
general direction, covering other input states, other cases 
of realistic noise, and investigating physical situations 
where our findings may be relevant for experiments, such 
as quantum magnet ometry [9j 10\ . 
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Appendix A: Coefficients for the non-unitary map 

The coefficients a, 6, c, and /, for the map defining 
the non-unitary evolution considered in the main text, 
are real and depend on cj, 7, and t. They are given by 



5 — lg-^t(l+acc+ay-az)7 _^ g^(aa3+ay)7^ 



/ = 



2^(a,-a^)^72-4c^2 
[a^j, - ayf ^'^ -Auj'^ 



(Al) 
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Appendix B: Details on the channel extension 
method 



As a result from the extension of the input state space, 
from to 2N particles, another upper bound on QFI can 
be obtained from eq. (15) of the main text, which does 
not involve any input state optimisation [16]: 

TN{t) < 4 min {iV \\aK{t)\\ + N{N - 1) WMt)^}. 

(Bl) 

The ll'll denotes the operator norm and 



aK{t) = Y^kl{t)k,{t), 

i 

PKit) = iJ2kKt)K,it) 



(B2) 



with Ki{t) = du;Ki{t). The minimization in (Bl) is per- 



formed over locally equivalent Kraus representations [15] 
of These are generated by all rxr Hermitian matri- 
ces h(t), so that K{t) = K{t) and K(t) =K(t)-ih(t)K(t), 
where K(t) is a column vector containing any starting 
{Ki{t)Y-^-^ as its elements and [hK.]i = '^j hijKj. Con- 
structing matrices 



Mt) = 

m = 



VKt2 K(t)t 
VX~a2 (iK(t)tK(t)) 

iK(t)tK(t) y/X~bl2 



(B3) 



with Id representing a d x d identity matrix, and Xa/b 
bei ng s ome real parameters, the minimisation problem 
of (Bl) is equivalent to a semi-definite programming 



task. This is achieved by realizing that the positive semi- 
definiteness of A{t) and B(t) corresponds respectively to 
the conditions 



aK{t) = K\t)K{t)<Xat2. 
f^K{t)f^K{t) = K\t)K{t)K\t)K{t) < A5I2, 



(B4) 



so that ( |B1[ ) can be rewritten as in eq. (16) of the main 
text. This construction has been introduced in [23] as the 
finite- A/" CE method, proving that an upper bound on J^n 
can be efficiently computed for any N. Its asymptotic 
version leading to a weaker bound has been described 
in and corresponds to adding a further constraint on 
h(t), namely that /3x(^) = 0, so that A5 = and the second 
term in eq. (16) vanishes. However, as highlighted in the 
main text, in our analysis we need to be more careful, 
as the ^-optimisation results in t being a function of N. 
For fixed t the asymptotic method of Ref. [15 directly 
applies, so that by making the constraint on h(t) that 
Pxit) = 0, from the first term of bound (Bl), we obtain 



the asymptotic constant factor SQL improvements for 
our model. 

For the parallel dephasing, as the asymptotic CE 
method coincides with the geometric method of fig. 2 [15], 
we obtain the weaker bound of eq. (6). In the case of the 
transversal noise, based on numerical results, we deter- 
mine the short-times asymptotic SQL constant factor to 
be Cxij^t) = ^t'^ -\-0{t^) what improves the analytically 
computed geometrical bound of eq. (8). However, if one 
substitutes in this formula the topt ( A^) evaluated with the 
finite-N CE method (stated above eq. (10)), one would 
arrive at a non-saturable asymptotic bound below the 
correct one-(lO) by a factor of 3. This fact emphasizes 
the applicability of finite-N CE method, which should al- 
ways be used to construct best asymptotic bounds when 
performing t-optimisation and considering channels of 
asymptotic SQL constant factor unbounded from below 
as t^O. 



